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The mathematics of wave propagation in a viscoelastic filament subjected to constant velocity 
transverse impact is discussed. The equations governing the stress-strain-time behavior are as- 
sumed to be those for a linear model consisting of a spring coupled in parallel with a spring and dash- 
pot in series. The nature of the solution is discussed, and a method is described for calculating the 
configuration, and the stress, strain, and particle velocity distributions along the impacted filament. 
The method used consists of an integration along the characteristics of the system of differential 
equations describing the problem. 
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1. Introduction 

When a flexible filament is struck transversely by 
a high velocity projectile, strain waves and waves of 
transverse motion are generated. The associated 
displacement of filament material can be recorded 
by high speed photography, and the resulting data 
can be analyzed to obtain the stress-strain, creep, and 
stress relaxation behavior of the material under im- 
pact conditions. Research of this type has been 
described by Petterson et al., [1, 2, 3] 1 and Smith 
etal.,[4,5,6]. 

A theory originated by Rachmatulin [7], Taylor 
[8, 9], and von Karman [10], and modified by Smith 
et al., [11, 12] has been used to analyze the experi- 
mental data. This theory, however, requires addi- 
tional modification in order to account adequately 
for viscoelastic effects. Malvern [13, 14] has helped 
in this respect by proposing a simple rate-dependent 
constitutive equation for approximating these effects, 
and has used it in treating the longitudinal impact 
problem. Further insight may be gained through the 
work of Lee and Kanter [15], Glauz and Lee [16], 
Morrison [17], and Smith [18] who have investigated 
the behavior of various linear spring and dashpot 
models subjected to longitudinal impact. 

Transverse impact of materials having rate-de- 
pendent constitutive equations has received less 
attention. Cristescu [19, 20, 21, 22, 23] has discussed 



1 Figures in brackets indicate the literature references at the end of this paper. 



various aspects of the problem of strain wave propa- 
gation and transverse wave motion in strings, including 
some strings with stress-strain-time behavior given 
by rate dependent constitutive equations. For some 
cases solutions calculated by a digital computer were 
obtained by integration along the characteristics of 
the system of differential equations describing the 
problem, but usually only continuous solutions were 
considered. 

The problem of transverse impact at constant 
velocity involves a solution in which the strain is dis- 
continuous at the strain wave front, and the slope of 
the filament is discontinuous at the front of a wave of 
transverse motion. This latter problem will be dis- 
cussed here, and a solution obtained by an extension 
of the methods of Cristescu. The understanding 
gained is expected to be of use in the design and 
analysis of future experiments. 



2. The Model 

The constitutive equation governing the stress- 
strain-time behavior of the filament is that for the 
three-element linear spring and dashpot model 
depicted in figure 1. The longitudinal impact behavior 
of this model has been treated previously by Morrison 
[17] and Smith [18]. In the model the Hooke's law 
constants of the two springs are g\ and £2. and the 
viscosity constant of the dashpot j\. T is the total 
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Figure 1. 



Linear spring and dashpot model consisting of two 
springs and one dashpot. 



nates respectively. At time £ = let the filament at 
x —Q be impacted at constant velocity V in the y 
direction. The motion of the filament is to be de- 
scribed in terms of a Lagrangian coordinate system; 
i.e., a system of x coordinates fixed to the filament, 
moving and extending with it, but such that in the 
initial unstrained state the x and x coordinates of the 
filament coincide. The displacement of any point x 
of the filament at time t relative to its original position 
is given by £(x, t) for horizontal displacements, and 
7)(x, t) for vertical displacements. The location in 
the laboratory coordinate system of a point x on the 
filament is thus given by 



x' =x-\-Jj(Xi t), y' = rj{x,t). 



(3a, b) 



Consider an element of unstrained length A*. After 
impact, this element is strained, and its new length 
measured in the laboratory coordinate system becomes 

Ax V( l + d£/d x) 2 -\- (dr)ldx)' 2 . The increase in 



per unit length, or strain e, is thus 



6= VTT-f di;ldx) 2 + (d^/dx)' 1 - 1 . 



length 



(4) 



The angle 6 that the element makes with the hori- 
zontal is given by 



force, e the strain, and t the time, 
equation is 



The constitutive 



i I \ ^€ glg 2 BT gl rp 

dt yi dt y, 



(1) 



It is more convenient to represent the model in 
terms of other parameters. Thus let gi = kg, g2 = 
(1 — \)g, y\ = \gT. The parameter r = y\lg\ is the re- 
laxation time of the Maxwell element constituting one 
branch of the model, and X is a dimensionless parame- 
ter that takes on values between and 1. When 
A=0 the model degenerates into a spring of modulus 
g, and when k=\ the model becomes a Maxwell 
element. For practical applications a value of \ of 
approximately 0.1 coupled with a relaxation time r 
less than the observation time following an impact is 
of interest. Such a model would be mostly springlike 
in its behavior, but would exhibit noticeable visco- 
elastic effects shortly after impact. Another con- 
venient parameter is c= Vg/m, the velocity at which a 
strain wave front propagates in the filament. The 
quantity m is the mass per unit length of the unstrained 
filament. Equation (1) expressed in the parameters 
r, \, c becomes 



de , (1 —\)mc 2 

1 6 : 

dt r 



dt T 



(2) 



3. Formulation of the Problem 

In the laboratory or observer's coordinate system, 
let x and y' denote the horizontal and vertical coordi- 



cos z 



1 



1 + e 



1 + 



dx 



sin 0- 



1 dt) 
1+e dx 



(5a, b) 



so that the differential equations that govern the motion 
of the small element are 






dt- ' 



dx dx 



1+e 



1 + 



dl 
dx 



_d_ 
dx 



T sin 0) = — 
dx 



T dr) 



1+edx 



(6a, b) 



Equations (6a) and (6b) can be incorporated into a 
system of linear equations by introducing the hori- 
zontal and vertical components of the velocity (in x\ 
y space) of a point on the filament 



-dtldt. 



v = dr)ldt 



and the direction variables a, (3 defined by 

_d_ri 



a= \ H , 

dx 



P- 



dx 



(7a, b) 



(8a, b) 



It is also helpful to introduce the flow or particle 
velocity w of a point along the direction of the filament. 
An increment of this velocity dw has horizontal and 
vertical components du and dv respectively and thus 
satisfies the relation 



dw—du cos 6-\-dv sin 6. 



(9) 
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The system of linear equations becomes 

aT 



du _ d 
dt dx 

dv = d_ 
dt~ dx 

du _da 
dx~~ dt 



m(\+e) 



BT 



/M(l+€). 



(10a) 



(]()!)) 



(10e) 




dv = dfi 
dx dt 



*-ITI d " + i+7* 



dw = du + 7^7 — dv 

1 +e 1 +e 



(10d) 



(10e) 



(10f) 



FIGURE 2. Configuration of a spring-like filament at time t after 
transverse impact. 

The filament, originally horizontal, was struck with velocity V in the vertical direction. 
After time / the resulting strain wave propagating with velocity r has a front at F, and the 
filament becomes uniformly strained at strain €„. The strained material between points 
F and B moves inward horizontally with velocity —r« » but at the transverse wave front, 
point B, this inward motion ceases. Material in the transverse wave moves vertically 
with velocity V . The transverse wave front moves with velocity P n relative to points fixed 
on the filament, but with velocity (1 + e )/'<> — ee<» with respect to the laboratory (x\ y') 
coordinate system. 



mc 

dt 



(1 -k)mc 2 



dT 

dt 



T. 



(lOg) 



Equation (lOe) is the differential form of eq (4). 

The seven variables u. v, ilk a, /3, e, and T must 
satisfy the initial conditions: For t = 0, x > 

u = v=w= l3 = 6=T={). a=] (] ]) 

and the boundary conditions: For t > 0, x = () 

u = 0,v=V (12a, b) 

and for t > 0, #-*<» 



u = v- 



r=o, a=l. 



(13) 



If conditions (12) are substituted into eq (lOf), 
additional boundary condition: For t> 0, x = 



w= a constant 



(14) 



is obtained. The value of this constant velocity is 
attained at the instant of impact, before any creep or 
stress relaxation has occurred, and thus is the same as 
that given in the solution for the springlike filament 
discussed in the next section. 

In addition to the initial and boundary conditions, 
certain compatibility conditions, to be given later, 
must be satisfied. 

4. Solution for the Case A = 

When X=0 the model degenerates to a spring of 
Hooke's law constant g. The solution for this case, 
discussed in detail elsewhere [11] is illustrated in 



figure 2. The impact causes a wave of constant strain 
e and constant tension T {) to propagate along the 
filament at velocity c, so that at time t the strain wave 
front has arrived at point F. In the wake of this 
wave front the strained material flows inward at ve- 
locity — ceo. The increase in length of the filament 
is taken up by a wave of transverse motion, which at 
the time t has a front at point B. 

In the transverse wave the filament material moves 
in the y direction at the impact velocity V, and there 
is no component of velocity in the x direction. The 
profile of the wave is a straight line. The transverse 

wave front moves at velocity Po= V To/m{l + €o) rela- 
tive to the Lagragian coordinate system, so that the 
length of strained filament AB as measured in the 
x\ y' system is equal to (1 + €o)PqT, and the horizon- 
tal projection of this length is [(l + eo)Po — C€o]t' 

As OAB is a right triangle, the following relationship 
applies 

P = (H-eo) 2 Po 2 -[(l + e )Po-c€o] 2 . (15) 

If the mass per unit length m and Hooke's law con- 
stant g are known lor the specimen, the velocities 
c and Po can be calculated for any strain e , and upon 
substituting these values into eq (15). the value of 
€o corresponding to any transverse impact velocity 
V can be found. 

In the region between the strain wave front and the 
transverse wave front, the particle velocity w and the 
horizontal velocity u are both constant and equal to 
— ce<). In the region of the transverse wave. u = but 
the velocity w remains equal to — ceo, even though in 
the x\ y coordinate system there is no observable 
flow of material longitudinally along the filament. 
This paradox can be explained by the following inter- 
pretation of the particle velocity w. Suppose that 
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there is a second filament, inextensible and fixed to 
the impacted filament at a point a great distance in 
advance of the strain wave front. Before impact let 
the two filaments coincide with the inextensible fila- 
ment ending at the point of impact, and after impact 
let the inextensible filament be constrained to move 
alongside the impacted filament. Under these con- 
ditions the point under impact will recede from the 
end of the inextensible filament with velocity w. 

The solution for the limit case \=0 can be stated 
briefly as follows. The quantities w, e, and T are 
continuous between the point of impact and the strain 
wave front at x = ct, and have the following values: 



w — — ceo, € — €o, T= 



(16a,b,c) 



where € is found by solving eq (15). The quantities 
w, i;, a, /3 are discontinuous at the transverse wave 



front x = P t = A/r /77z(l+eo). 
For ^ x < P t 

u = 0, v=V 



These equations then can be arranged to obtain 

(dw dw\ (dT dT\ 1 T 21 

\dt dx J \dt dx) r 

(dw dw\ (dT dT\ 1 2 



For differentiations along curves of the families 
given by 



dx/dt z 



(20) 



the expressions in parentheses may be replaced by 
total derivatives so that the equations reduce to: Along 
the curves given by dx/dt = c 

mcdw -dT=[T-(l- \)mc 2 e ] dt/r (2 1 a) 

and along the curves given by dx/dt = — c 

mcdw + dT=- [T-(\-\)mc 2 e]dt/T. (21b) 



a=(l+€o) cos = H-€o — C€o/Po(17a,b,c,d) The curves satisfying eq (20) are known as charac- 
teristics. 
0=(l + €o) sin = -V/P o If /'= V77m(l + € ), eq (10a) can be expressed in the 

i r r> form 

and lor r t ^ x < ct 



u = w = — c€o* v=P= 0. a= 1 +€o 



(18a,b,c) 



The limit solution just presented indicates the type 
of solution to be sought for the general case where 
X > 0: namely, a solution in which w* e, and T are 
continuous between the point of impact and the strain 
wave front, and u, v, a, and (3 are continuous only be- 
tween the point of impact and the transverse wave 
front. A method for finding this solution by integrating 
along the characteristics of the system of differential 
equations will be developed in the succeeding sections. 

5. The Characteristics of the System of 
Equations 

If the variables u, v. a, and (3 are eliminated from 
eqs (10a) through (lOg) the following set of equations 
results. 2 





dw_ 1 dT 






dt m dx 






dw _de 
dx ~ dt 






i d€ (1- 

mc + 
dt 


-k)mc 2 

r 



e = f+ l -T. 

dt r 



(19a,b,c) 



2 Equations (19a) through (19c) plus the boundary conditions of eqs (11) through (14) 
involving w. e. and T are the same as the equations and boundary conditions for longitudinal 
impact at constant velocity —ceo. The distribution of w. e. and T can be calculated from 
the solutions for this case given by Morrison [17| and Smith [18]. However, the subsequent 
determination of values for //. r. a, and /3 from eqs (10a) through (lOd) would be difficult 
and will not be attempted here by this method. 



dt dx dx 



which can be combined with eq (10c) to obtain the 
expressions 

dt dx \dt dx * dx 



du du 

dt dx 



\dt dx / dx 



Along curves of the families given by 
dx/dt = ±P 



these expressions reduce to 

dt dt dx 



Similarly eqs (10b) and (lOd) yield 

dt dt H dx J 



(22) 



(23a) 



(23b) 



Equations (23a) and (23b) can be combined to give: 
Along the curves given by dx/dt = P 



fidu- adv = Ppda- Padfi 



(24a) 
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tities with respect to x and /, however, are not nec- 
essarily continuous at the strain wave front and at 
(24b) tne transverse wave front. At these places of dis- 
continuity, according to the principle of kinematic 

The curves satisfying eq (22) are additional charac- compatibility [24] , the following relations hold: 

teristics. 



and along the curves given by dx/dt = — P 
fidu -adv = - Ppda + Pad/3 



Other equations holding along these characteristics 
can be found by the following manipulation involving 
eqs (19a) and (19b). 

d_ (1 +6) (dTldx)-T(deldx) 

dx { ) m(l+e) 2 



dx 

drj 

dx 



dx 
dt* 



dx 
dt 



+ 



dt 

~dt 



= 



(27a,b) 



1-he 



and 



]_dT 
m dx 



dw 
Yt 



dw 
dt 



,J ™ = ^ _ p 



dx 



1 dT 

m dx 



dx 



But along curves defined by dx/dt = P 



dx dt dt 



fJ (k_ p dw 



dt 



dx 



where the bracketed quantities represent jumps in 
value occasioned by the passage of a wave front at 
the point under consideration. For instance, at the 
point x the expression [e] is denned as 



[e] = lim {e(x, t + d)- 

8^0 



€(X, t)}. 



The strain wave front velocity is greater than the 
transverse wave front velocity. Therefore at the 
strain wave front, where dx/dt = c, it follows that 17 = 
and, from eq (4), d^/dx = e. Equation (27a) then yields 
the relation 



that 



dx 



1-fe 



dw p de 

lt~ dt 



Thus along curves defined by dx/dt = P, eqs (23a) and 
(23b) reduce to 



( 1 + e)(du - Pda)= a( dw - Pde) 
( 1 + e)(dv - Pdfi) = (3( dw - Pde ) 

Similarly along curves defined by dx/dt = — P 
(1 + e)(du + Pda)=a(dw + Pde) 
( 1 + e)(dv + Pdfi) = /3(dw + Pde) 



(25a,b) 



(26a,b) 



Along characteristic curves in x, t space the inte- 
gral surface of the system of differential equations 
may have discontinuities. It should be noted that the 
characteristic curve x = ct determines the discon- 
tinuous wave front for the stress, strain and particle 
velocity distributions, arid the characteristic curve 
that passes through the origin and is a solution of the 
equation dx/dt = P, determines the front of the trans- 
verse wave where the quantities u, v, a, and f3 are 
discontinuous. 

6. Compatibility Conditions 

An obvious requirement of the solution is that the 
displacements £ and 77 be continuous at all points along 
the filament. The partial derivatives of these quan- 



[e]=-- c [w]. 



(28) 



Similarly at the transverse wave front, where dx/dt = P, 
eqs (27a) and (27b) yield 



[«-l]=~[u] 



L3]=- -plv] 



(29a,b) 



As [a— 1] =a— 1 — e, [u]=u — w, [fi] = P, and 
[u]=i>, where a, u, /3, and v refer to values in the 
wake of the wave front, eqs (29a) and (29b) can be 
expressed more conveniently as 



a=l + €-ulP + wlP 
p=-v/P 



(30a,b) 



Another compatibility condition is obtained as fol- 
lows. Consider a point x on the filament. Just before 
arrival of the strain wave front at time t, the tension 
T and particle velocity w are both zero. At a slightly 
later time t-\- At _=t~\- Ax/c, the tension is — Ax(dT/dx) 
and the particle velocity is (Ax/c) (dw/dt). Thus the 
jump relations for T and w are \T\ = — Ax(dT/dx) and 
: (Ax/c)(dw/dt). Application of eq (19a) then 



w 



yields the condition of dynamic compatibility 
[T\ = —mc[w\. 



(31) 
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7. Values of the Tension, Strain, and Particle 
Velocity at the Strain Wave Front 

Equation (21a) is a differential relation involving 
the tension T, strain €, and particle velocity w, that 
holds at the strain wave front. Hence if values of T 
and w from eqs (28) and (31) are substituted into this 
relation, it becomes 2de=—{\elr)dt, which when 
integrated yields 



e(ct, t) = e(0,0)e-( kt ' 2T \ 



(32) J2 

X 

< 



Substituting this result into eqs (28) and (31) then 
gives 



w(ct, t) = -ce(0,0)e- {kt ^ 
T(ct, t) = mc 2 e(0,0)e-WT ) . 



(33) 
(34) 



When x = t = the strain and transverse wave fronts 
coincide. If values at t = x = are designated by a 
zero subscript, the boundary conditions give ^o = 0, 
vo — V, and eqs (30a) and (30b) reduce to 




x AXIS 



FIGURE 3. A plot of the characteristics network. 

The straight line characteristics of positive slope are solutions of the equation dx/dt = c. 
The straight line characteristics of negative slope are solutions of the equation dx/dt = — c. 
The straight line characteristic through the origin has the equation x = ct, or equation of 
the strain wave front. The curved ch aracteristic passing through the origin is that solution 
of the equation dx/dt= VT/mQ +e) that gives the location of the transverse wave front. 



a = 1 + e — ceolPo 
Po=-V/Po 



(35a,b) 



Thus eq (4) expressed as al + fil = (l + e ) 2 , can be put 
in the form V 2 = P§(1 + e ) 2 - [(1 + € )P - C€ ] 2 . But 
Po = c Ve /(l + €o) so that 



P/c 2 = eo(l + e )-[Veo(l + eo)-€o] 2 . 



(36) 



The value of e to be used in eqs (32), (33), and (34) 
is obtained by solving eq (36). 

8. Integration Along the Characteristics 

A numerical solution to the system of differential 
eqs (10a) through (lOg) and the initial and boundary 
conditions (11) through (14) can be found by integrating 
along the characteristics of the system. The nature 
of the problem and the method of solution are easily 
understood from an inspection of figure 3, in which 
some of the characteristics are plotted. 

The straight line characteristics of positive slope are 
solutions of the equation dxjdt = c. Their slope is 
equal to 1/c. The other straight line characteristics 
completing the network are solutions of the equation 
dxjdt = — c. The characteristic passing through the 
origin, and having the equation x — ct corresponds to 
the discontinuity marking the wave front of the strain, 
tension, and particle velocity distributions. Thus in 
the domain between this characteristic and the x axis 
the strain e, tension T, and particle velocity w are zero. 
Along the characteristic the values of e, w, and T are 
given by eqs (32), (33), and (34). These equations thus 



may be thought of as boundary conditions holding 
along the characteristic x = ct. 

The curved line passing through the origin represents 
a solution of the characteristic equation dx/dt = P. 
The exact course of this characteristic is not known in 
advance, as P is a function of e and T. This charac- 
teristic corresponds to the front of the wave of trans- 
verse motion that propagates along the filament; thus 
in the domain between this characteristic and the 
characteristic x = ct the velocity u is equivalent to the 
velocity w, the direction variable a is equal to 1 + e, 
and the velocity v and direction variable /3 are equal to 
zero. The jump relations (30a) and (30b) holding along 
the transverse wave front characteristic constitute a 
second set of boundary conditions. 

A third set of boundary conditions expressed by eqs 
(12a), (12b), and (14) holds along the time axis. 

The problem consists in finding the values of w, 
e, and T at all of the mesh points in the domain be- 
tween and along the time axis and the characteristic 
x = ct, in finding the course of the transverse wave 
front characteristic, and in finding the values of u, 
v, a, and f3 at all of the mesh points between and 
along the time axis and the transverse wave front 
characteristic. The values of these quantities must 
satisfy the differential relations (lOe), and (lOf), the 
relation (lOg) holding along the line x = a constant, 
the differential relations (21a), (21b), (24a), (24b), 
(25a), (25b), (26a), and (26b) holding along the char- 
acteristics, and the boundary conditions. 

A first step towards the solution consists in solving 
eq (36) for the initial strain e at the point of impact. 
The values of e, w, and T along the characteristic 
x = ct are then given by eqs (32), (33), and (34). 

Values of e, w, and T are then found at the first mesh 
point on the time axis above the origin. Figure 4a 



262 




Figure 4. Mesh situations encountered during integration of the 
system of differential equations along the characteristics network. 



depicts one of the triangular meshes lying along the 
t axis. In the lowest mesh points the values of 
e, w, and T are known, and from the boundary condi- 
tion (14) w is known to be constant along the t axis. 
The strain 6 and tension T are found by solving eq 
(21b), which holds along the characteristic of negative 
slope, and eq (10g), which holds along the t axis. 
These equations can be expressed in the form 

T- T> = - [f 2 -(1 - k)mc 2 E 2 ]&tlT - mdwx - w*) 

mc 2 (e-€i)-(T-T 1 ) = 2[T 1 -(l-\)mc 2 E 1 ]AtlT 



(37a, b) 



and 4b, values of e, iv, and T can be found at each <>i 
the mesh points in the" domain between the / axis 
and the characteristic x = ct. 

The values of a and /3 at the origin are next calcu- 
lated from eqs (30a) and (30b). The values of // and v 
are already known from the boundary conditions. 
The solution for the mesh situation depicted in figure 
4c can then be used to determine the course of the 
transverse wave front characteristic and the values of 
e , w, 7\ u, v, a , and (3 at its intersections with the other 
characteristics. In figure 4c the quantities €, w, T, u, 
v, a, and /3 are known at point 1, and the values of 
e, w, and T are known at points 2 and 3. The values 
of e, w* 7\ u, v, «, and /3 are required at the intersec- 
tion of the transverse wave front characteristic with 
the characteristic passing through points 2 and 3. 
As a first step P is set equal to v Ti/m(l -h €i) and the 
line (x—Xi) = P(t — ti) is drawn through point 1 to 
intersect the characteristic between points 2 and 3. 
The values of T and e at the intersection are deter- 
mined by interpolation. The quantity P is then set 

equal to vTlm{\ +E), where T and E are the averages 
of T and 7\, and € and ei respectively. One or two 
iterations are performed to obtain more accurate values 
of T and e and the coordinates (x 9 i) of the intersection. 
The value of w at the intersection is then found by 
interpolation. The values of u, v, a, and /3 at the 
intersection are found by solving the set of four equa- 
tions (30a), (30b), (25a), and (25b) expressed in the 
form 

A 

u — u\ — P(ol — cxi)= , . r [w — wi — P{e — €i)\ 

1 -+- 1 

u + Pa = w + P(l+e) 



where At=t — t2 and 2At = t — t\. An iter ative_ proc- 
ess is used in the solution. The quantities T\, T 2 , 
E\, and E 2 are first set equal to T\, T>>, €i, and e->, re- 
spectively, and approximate _yalues are found for T 
and e. The quantities Ti, T>>, E u and E 2 are then 
assigned the average values (7 7 + Ti)/2, (T-\-T->)/2. 
(€ + €i)/2, and (e + e 2 )/2, respectively, and eqs (37a) 
and (37b) are solved again for more accurate values 
of T and e. Another iteration may be performed if 
desired. 

Figure 4b depicts a rhombus shaped mesh formed 
by the intersection of characteristics. If values of 
e. w. and T are known at each of points 1. 2, and 3 
the values of e, w, and T at the vertex are found by 
solving the system of equations; eq (21a) which 
holds along the upper characteristic passing through 
point 1, eq (lOg) which holds along the vertical line 
passing through point 2. and eq (21b) which holds 
along the upper characteristic passing through point 3. 
The solution is found by an iterative process similar to 
that just described. 

By applying successively the methods of solution 
associated with the mesh types depicted in figures 4a 



v + P[3 = 0. 



(38a,b,c,d) 



In the solution process A and B are first set equal to 
ol\ and ft i respectively and in subsequent iterations 
set equal to the average values along the transverse 
wave front characteristic. 

Solutions are needed for two additional mesh situa- 
tions in order to find the values of u, v, a, and (3 at the 
mesh points between the t axis and the transverse 
wave front characteristic. Figure 4d depicts a mesh 
of the type found along the t axis. The characteristic 
between points 1 and 3 may be either the transverse 
wave front characteristic or one of the solutions of the 
characteristic equation dx/dt = c. Usually the mesh 
is triangular in shape with point 1 lying on the t axis, 
but sometimes when the transverse wave front charac- 
teristic has a steep slope, the mesh has the more 
general form shown. The values of e, w, T, u, v, a, /3 
and the coordinates x and t are known at points 1 and 3, 
and the values of u, v 9 w, e, and T are known at point 4. 
The values of a and (3 at point 4 are required. As a 
first step the characteristic of slope — 1/P is drawn 
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between points 4 and 2. In this construction iteration 
and interpolation processes are used to determine 
values of P, e 2 , and T 2 , and the coordinates (x 2 , t 2 ) of 
point 2. Interpolated values for u 2 , i>2, a 2 , 02, and w 2 
are calculated next. The values of a and /3 at point 
4 are then obtained by solving eqs (26a) and (26b). 
Figure 4e depicts a four-sided mesh in which the 
characteristic between points 3 and 5 may be either 
the transverse wave front characteristic or a solution 
of the characteristic equation dx/dt = c. Values of 
u, v, w, a, /3, e, and T are known at points 1, 3, and 5, 
and values of e, w, and T are known at point 6. The 
values of u, v, a, and j8 at point 6 are required. To 
find the solution lines of slope + 1/P and -1/P are 
constructed to intersect at points 2 and 4, respec- 
tively, and values of u, v, w, a, p, e, and T are deter- 
mined at these points. The values of u and a at point 
6 are then obtained by solving eqs (25a) and (26a) 
expressed in the form 

u — uo — P(a — a>>) = ; , r [w s — w 2 — P(e 6 — € 2 )] 

1 TL 



ua + P(a _ ola) '- 



l+E 



[w G — W4 + P(t6 ~ €4)] • 

(39a,b) 



The values of t; and /3 at point 6 are obtained similarly 
by solving eqs (25b) and (26b). 

By systematically applying the methods of solution 
outlined for the five mesh types depicted in figure 4, 
it is possible to calculate numerical values of e, w, 7\ 
a, v, a, and (3 at any distance x along the filament, 
and at any time t. 

In order to find the configuration of the filament in 
laboratory (x r , y) coordinates it is only necessary to 
calculate appropriate values of f (x, t) , and 77 (x, t) by 
integrating u{x, t) and v(x, t) with respect to time. It 
is of special interest to perform this calculation in 
order to determine how much curvature is introduced 
into the transverse wave portion of the filament by the 
creep and stress relaxation processes that occur after 
impact. 



The methods of calculation just described have been 
incorporated into a set of programs for use with a digital 
computer. Calculations have been made for various 
values of the propagation velocity c, impact velocity V, 
relaxation time r, and parameter X that approximate 
conditions realized in tests on real materials. It is 
planned to present some of the results of these calcula- 
tions in a later publication. 
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